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ABSTRACT 

The subject of this paper is a special class of three-degree- 
of-freedom parallel manipulators. The singular configurations of 
the two Jacobian matrices are first studied. The isotropic con- 
figurations are then found based on the characteristic length of 
this manipulator. The isoconditioning loci for the Jacobian ma- 
trices are plotted to define a global performance index allowing 
the comparison of the different working modes. The index thus 
resulting is compared with the Cartesian workspace surface and 
the average of the condition number. 



I Introduction 

Various performance indices have been devised to assess 
the kinetostatic performances of serial and parallel manipulators. 
The literature on performance indices is extremely rich to fit in 
the limits of this paper, the interested reader being invited to look 
at it in the references cited here. A dimensionless quality index 
was recently introduced in [1] based on the ratio of the Jaco- 
bian determinant to its maximum absolute value, as applicable to 
parallel manipulators. This index does not take into account the 
location of the operation point of the end-effector, because the 
Jacobian determinant is independent of this location. The proof 
of the foregoing result is available in [2], as pertaining to serial 
manipulators, its extension to their parallel counterparts being 
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straightforward. The condition number of a given matrix, on the 
other hand, is well known to provide a measure of invertibility of 
the matrix [3 ]. It is thus natural that this concept found its way in 
this context. Indeed, the condition number of the Jacobian ma- 
trix was proposed in El as a figure of merit to minimize when 
designing manipulators for maximum accuracy. In fact, the con- 
dition number gives, for a square matrix, a measure of the relative 
roundoff-error amplification of the computed results 1 3 ] with re- 
spect to the data roundoff error. As is well known, however, the 
dimensional inhomogeneity of the entries of the Jacobian matrix 
prevents the straightforward application of the condition number 
as a measure of Jacobian invertibility. The characteristic length 
was introduced in [5 ] to cope with the above-mentioned inhomo- 
geneity. 

In this paper we use the characteristic length to normalize 
the Jacobian matrix of a three-dof planar manipulator and to cal- 
culate the isoconditioning loci for all its working modes. 



2 Preliminaries 

A planar three-dof manipulator with three parallel PRR 
chains, the object of this paper, is shown in Fig. 1. This ma- 
nipulator has been frequently studied, in particular in (6]|3 • The 
actuated joint variables are the displacements of the three pris- 
matic joints, the Cartesian variables being the position vector p 
of the operation point P and the orientation of the platform. 

The trajectories of the points A, define an equilateral triangle 
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Figure 1 . A three-DOF parallel manipulator 



whose geometric center is the point O, while the points B\, # 2 
and #3, whose geometric center is the point P, lie at the corners 
of an equilateral triangle. We thus have a/ = % + (z — l)(27i/3), 
for i= 1,2,3. Moreover, / = l\ = Z 2 = h, with /; denoting the 
length of AiBi and r = r\ = r 2 — r 3 , with n denoting the length 
of BiP, in units of length that need not be specified in the paper. 
The layout of the trajectories of points A; is defined by the radius 
R of the circle inscribed in the associated triangle. 

2.1 Kinematic Relations 

The velocity p of point P can be obtained in three different 
forms, depending on which leg is traversed, namely, 



p = ai+f|iE(bi -ai) + 6E(p-bi) 
p = a 2 + f|2E(b 2 - a 2 ) + 0E(p - b 2 ) 
p = a 3 + f| 3 E(b 3 - a 3 ) + 6E(p - b 3 ) 



(la) 
(lb) 
(lc) 



with matrix E defined as 



E 



The velocity a ; of points A\ is given by 
^ _ Pl IIP.II — p * L sLn(«,) 



p/e,- 



We would like to eliminate the three idle joint rates f|i, f| 2 
and f| 3 from eqs.([T^-c), which we do upon dot-multiplying the 
former by (b; — a/) 7 , thus obtaining 



(bi - ai) Y p = (bi - ai) Y piei + (bi - ai) Y 0E(p - bi) (2a) 

(b 2 - a 2 ) r p = (b 2 - a 2 ) r p 2 e 2 + (b 2 - a 2 ) r eE(p - b 2 ) (2b) 

(b 3 - a 3 ) r p = (b 3 - a 3 ) r p 3 e 3 + (b 3 - a 3 ) r 0E(p - b 3 ) (2c) 

Equations d2^-c) can now be cast in vector form, namely, 



At Bp with t : 



and p 



Pi 
P2 
P3 



(3) 



with p thus being the vector of actuated joint rates. 

Moreover, A and B are, respectively, the direct-kinematics 
and the inverse-kinematics matrices of the manipulator, defined 
as 



B 



(bi-ai) r 
(b 2 -a 2 ) r 
(b 3 -a 3 ) r 

(bi -ai) r ei 





(bi-ai) r E(p-bi) 
(b 2 -a 2 ) r E(p-b 2 ) 
(b 3 -a 3 ) r E(p-b 3 ) 


(b 2 -a 2 ) r e 2 

(b 3 -a 3 ) r e 3 



(4a) 



(4b) 



where e; is a unit vector in the direction of the /th prismatic joint. 



When A and B are nonsingular, we obtain the relations 

t = Jp, with J = A _1 B and p = Kt 
with K denoting the inverse of J. 

2.2 Parallel Singularities 

Parallel singularities occur when the determinant of matrix 
A vanishes ||8]|9l. At these configurations, it is possible to move 
locally the operation point P with the actuators locked, the struc- 
ture thus resulting cannot resist arbitrary forces, and control is 
lost. To avoid any performance deterioration, it is necessary to 
have a Cartesian workspace free of parallel singularities. For 
the planar manipulator studied, such configurations are reached 
whenever the axes A\B\, A 2 # 2 and A 3 # 3 intersect (possibly at 
infinity), as depicted in Fig. [2] 

In the presence of such configurations, moreover, the ma- 
nipulator cannot resist a force applied at the intersection point. 
These configurations are located inside the Cartesian workspace 
and form the boundaries of the joint workspace |8 |. 
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However, notice that B is dimensionally homogeneous, and 
need not be normalized. 



3.3 Isotropic Configuration 

In this section, we derive the isotropy condition on J to de- 
fine the geometric parameters of the manipulator. We shall obtain 
also the value L of the characteristic length. To simplify A and 
B, we use the notation 



Figure 2. Parallel singularity 



Figure 3. Serial singularity 



2.3 Serial Singularities 

Serial singularities occur when det(B) = 0. In the presence 
of theses singularities, there is a direction along which no Carte- 
sian velocity can be produced. Serial singularities define the 
boundary of the Cartesian workspace. For the topology under 



1, 


= (bi-a,- 


h 


= (bi-a,- 


nii 


= (bj-a,- 


Ji 


= LAiBiP 



(7a) 
(7b) 
(7c) 
(7d) 



We can thus write matrices A and B as 



study, the serial singularities occur whenever (b; — a i ) r e i = for 




[If 


-h/L 




m\ 





" 


at least one value of i, as depicted in Fig. [3] for i = 2. 


A = 




-k 2 /L 


B = 





mo 











-h/L 




_ 





m 3 _ 



(8) 



3 Isoconditioning Loci 

3.1 The Matrix Condition Number 

We derive below the loci of equal condition number of the 
matrices A, B and K. To do this, we first recall the definition 
of condition number of an m x n matrix M, with m < n, k(M). 
This number can be defined in various ways; for our purposes, 
we define k(M) as the ratio of the smallest, c s , to the largest, O/, 
singular values of M, namely, 



K(M) = J 



(5) 



The singular values {d",}™ of matrix M are defined, in turn, as 
the square roots of the nonnegative eigenvalues of the positive- 
definite mxm matrix MM 7 ". 



3.2 Non-Homogeneous Direct-Kinematics Matrix 

To render the matrix A homogeneous, as needed to define its 
condition number, each term of the third column of A is divided 
by the characteristic length L ifTUI . thereby deriving its normal- 
ized counterpart A: 



A = 



(bi- 
(b 2 - 
(b 3 - 



-*2) T 
-a 3 ) r 



(bi- 
(b 2 - 
(b 3 - 



-ai) r E(p-bi)/L 
-a 2 ) r E(p-b 2 )/L 
■a 3 ) r E(p-b 3 )/L 



(6) 



which is calculated so as to minimize k(A), along with the pos- 
ture variables pi, p2 and p3. 



Whenever matrix B is nonsingular, that is, when m\ ^ 0, for i 
1,2,3, we have 



K 



If /mi 
lf/m 2 

m/m 3 



-k\/{Lm\) 
-k 2 /(Lmi) 
-k$/(L m$) 



Matrix J, the normalized J, is isotropic if and only if KK 
x 2 13x3 for x > and K = J _1 , i.e., 



j/L')/mi = x' 

2 / T 2\ /, 2 _ ^2 



(l[l2 
(lfl3 
(lp3 



(l[li+^/L 2 )/m2 = x 2 

(l 2 l2 

{\lh+kl/L l )/mi = X 
+ kik 2 /L 2 )/(mim 2 ) =0 
+ k 1 k 3 /L 2 )/(mim 3 ) =0 
+ k 2 k 3 /L 2 )/(m 2 m 3 ) =0 



(9a) 
(9b) 
(9c) 
(9d) 
(9e) 
(9f) 



From eqs.(|9^-f), we can derive the conditions below: 



INI = INI = INI dOa) 

||p-bi|| = ||p-b 2 || = ||p-b 3 || (10b) 

l[l 2 = l£l 3 =l£l 3 (10c) 

m\in 2 = m\m 3 = m 2 m 3 (10d) 

In summary, the constraints defined in the eqs.JToh-df are: 
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o Pivots Bi should be placed at the vertices of an equilateral tri- 
angle; 

o Segments A{B[ form an equilateral triangle; 
o The trajectories of point A/ define an equilateral triangle an 
hence, / = l\ = h = h> 

Notice that the foregoing conditions, except the second one, were 
assumed in §2. 

3.4 The Characteristic Length 

The characteristic length is defined at the isotropic configu- 
ration. From eqs.([9]i-f), we determine the value of the character- 
istic length as, 




By applying the constraints defined in eqs.dg^-d), we can write 
the characteristic length L in terms of angle y, i.e. 

L = v // 2rsin(y) 

where y=J\ = J2 = 73, was defined in eq.(|7]i) and y G [0 2n] . 

This means that the manipulator under study admits several 
isotropic configurations, two of which are shown in Figs. 4a and 
b, whereas the characteristic length L of a manipulator is unique 
0. When y is equal to n/2, Fig. 4a, i.e., when A/5/ is perpendic- 
ular to BiP, the manipulator finds itself at a configuration furthest 
away from parallel singularities. To have an isotropic configura- 
tion furthest away from serial singularities, we have two condi- 
tions: ef E(b/ - a,-) = and r = R/2. 




Figure 4. Two isotropic configurations with two values of p 



3.5 Working Modes 

The manipulator under study has a diagonal inverse- 
kinematics matrix B, as shown in eq.(5b), the vanishing of one 
of its diagonal entries thus indicating the occurrence of a serial 
singularity. The set of manipulator postures free of this kind of 
singularity is termed a working mode. The different working 
modes are thus separated by a serial singularity, with a set of 
postures in different working modes corresponding to an inverse 
kinematics solution. 

The formal definition of the working mode is detailed in [0. 
For the manipulator at hand, there are eight working modes, as 
depicted in Fig. 5. 




(5) (6) (7) (8) 



Figure 5. The eight working modes of the 3PRR manipulator 



However, because of symmetries, we can restrict our study 
to only two working modes, if there are no joint limits. Indeed, 
working mode 1 is similar to working mode 5, because for the 
first one, the signs of the diagonal entries of B are all negative, 
and for the second are all positive. A similary reasoning is appli- 
cable to the working modes 2-6, 3-7 and 4-8; likewise, the work- 
ing modes 3-4 and 7-8 can be derived from the working modes 
2 and 6 by a rotation of 120° and 240°, respectively. Therefore, 
only the working modes 1 and 2 are studied. We label the corre- 
sponding matrices as A/, B/, K/ for the ith working mode. 

3.6 Isoconditioning Loci 

For each Jacobian matrix and for all the poses of the end- 
effector, we calculate the optimum conditioning according to the 
orientation of the end-effector. We can notice that for any orien- 
tation of the end-effector, there is a singular configuration. 

Figure depicts the isoconditioning loci of matrix A. 

We depict in Fig.[7]the isoconditioning loci of matrix B. We 
notice that the loci of both working modes are identical. This is 
due to both the absence of joint limits on the actuated joints and 
the symmetry of the manipulator. For one configuration, only 
the signs of m/ change from a working mode to another, but the 
condition number K is computed from the absolute values of ra/. 
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(a) (b) 



Figure 6. Isoconditioning loci of the matrix (a) Ai and (b) A2 with R/r = 

2 and l/r = 2 




Figure 7. Isoconditioning loci of the matrix (a) Bi and (b) B2 with R/r = 

2 and l/r = 2 



matrices A and K decreases regularly around the isotropic con- 
figuration. The isoconditioning loci resemble concentric circles. 
However, for the second working mode, the isoconditioning loci 
of both matrices A and K resemble ellipses. 

The characteristic length L depends on r. Two indices can 
be studied according to parameter R\ (i) the area of the Cartesian 
workspace, called 5, and (ii) the average of the conditioning, 
called k. 

The first index is identical for the two working modes. Fig- 
ure [9] depicts the variation of S as a function of R/r, for l/r = 2. 
Its maximum value is reached when R/r = 0.5. 

S * , 




-I — , — , — , — , — , — , — , — . — ^-L> 

0.5 1.5 2.5 3.5 4.5 R/r 

Figure 9. Variation of S as a function of R/r, with l/r = 2 



For the three matrices studied, K can be regarded as a global 
performance index. This index thus allows us to compare the 
working modes. Figsures fTQl [TT1 and fT2l depict k(A). k(B) and 
k(K), respectively, as a function of R/r, with l/r = 2. 



The shapes of the isoconditioning loci of K (Fig. 8) are sim- 
ilar to those of the isconditioning loci of A; only the numerical 
values change. 




-2.5 -2.0 -1.5 -1.0 -0.5 



0.0 0.5 

) 



1.0 1.5 2.0 2.5 x/r -2.5 -2.0-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 x /r 
(b) 



Figure 8. Isoconditioning loci of the matrix (a) Ki and (b) K2 with 
#/r = 2and//r = 2 



For the first working mode, the condition number of both 




K(Ai) 



-7^— K(A 2 ) 



-1 1 1 1 1 r 

0.5 1.5 2.5 



I — 1 — R/r 



3.5 4.5 

Figure 10. k(Ai) and k(A2) as a function of R/r, with l/r = 2 



The value of k(Ai) increases with R. At the opposite, the 
maximum value of k(A2) is reached when R/r = 2. For both the 
working modes studied, k(Bi) and k(B2) are identical for R/r 
fixed. For the first working mode, the minimum value of k(Ki ) 
and the maximum area of Cartesian workspace S occur at differ- 
ent values of R/r. This means that we must reach a compromise 
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Figure 1 1 . k(Bi) and k(B2) as a function o\R/r, with l/r — 2 




Figure 12. k(Ki) and k(K2) as a function oiR/r, with l/r = 2 



under these two indices. For the second working mode, there is 
an optimum of k(K2) close to the optimum of S , for R/r = 2. 



4 Conclusions 

We produced the isoconditioning loci of the Jacobian ma- 
trices of a three-PRR parallel manipulator. This concept being 
general, it can be applied to any three-dof planar parallel manip- 
ulator. To solve the problem of nonhomogeneity of the Jacobian 
matrix, we used the notion of characteristic length. This length 
was defined for the isotropic configuration of the manipulator. 
The isoconditioning curves thus obtained characterize, for ev- 
ery posture of the manipulator, the optimum conditioning for all 
possible orientations of the end-effector. This index is compared 
with the area of Cartesian workspace and the conditioning av- 
erage. The two optima being different, it is necessary to find 
another index to determine the optimum values. The results of 
this paper can be used to choose the working mode which is best 
suited to the task at hand or as a global performance index when 
we study the optimum design of this kind of manipulators. 
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